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Executive summary 

A bulk-flow analysis and computer program for prediction of the static load perfonnance 
and dynamic force coefficients of angled injection, orifice-compensated hydrostatic / 
hydrodynamic thrust bearings have been completed. The product of the research is an 
efficient computational tool for the design of high-speed thrust bearings for cryogenic 
fluid turbopumps. The study addresses the needs of a growing technology that requires of 
reliable fluid film bearings to provide the maximum operating life with optimum 
controllable rotordynamic characteristics at the lowest cost. 

The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a 
set of bulk-flow mass and momentum conservation and energy transport equations. Mass 
flow conservation and a simple model for momentum transport within the hydrostatic 
bearing recesses are also accounted for. The bulk-flow model includes flow turbulence 
with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing 
recesses and film lands. The cryogenic fluid properties are obtained from realistic 
thennophysical equations of state. Turbulent bulk-flow shear parameters are based on 
Hirs’ model with Moody's friction factor equations allowing a simple simulation for 
machined bearing surface roughness. A perturbation analysis leads to zeroth-order 
nonlinear equations governing the fluid flow for the thrust bearing operating at a static 
equilibrium position, and first-order linear equations describing the perturbed fluid flow 
for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth- 
order flow field equations renders the bearing flow rate, thrust load, drag torque and 
power dissipation. Solution to the first-order equations detennines the axial stiffness, 
damping and inertia force coefficients. The computational method uses well established 
algorithms and generic subprograms available from prior developments. 

The Fortran90 computer program hydrothrust runs on a Windows 95/NT personal 
computer. The program, help files and examples are licensed by Texas A&M University 
Technology License Office. 

The study of the static and dynamic performance of two hydrostatic/hydrodynamic 
bearings demonstrates the importance of centrifugal and advection fluid inertia effects for 
operation at high rotational speeds. The first example considers a conceptual hydrostatic 
thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The 
large axial stiffness and damping coefficients of the bearing should provide accurate 
control and axial positioning of the turbopump and also allow for unshrouded impellers, 
therefore increasing the overall pump efficiency. The second bearing uses a refrigerant 
R134a, and its application in oil- free air conditioning compressors is of great 
technological importance and commercial value. 

The computed predictions reveal that the LTL bearing load capacity and flow rate 
increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial 



stiffness has a maximum for a recess pressure ratio -0.55, while the axial damping 
coefficient decreases as the recess pressure ratio increases. The computer results from 
three flow models are compared. These models are a) inertialess, b) fluid inertia at recess 
edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia 
model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the 
other hand renders the largest damping coefficients and inertia coefficients. The most 
important findings are related to the reduction of the outflow through the inner radius and 
the appearance of subambient pressures. 

The performance of the refrigerant hybrid thrust bearing is evaluated at two operating 
speeds and pressure drops. The computed results are presented in dimensionless form to 
evidence consistent trends in the bearing perfonnance characteristics. As the applied axial 
load increases, the bearing film thickness and flow rate decrease while the recess pressure 
increases. The axial stiffness coefficient shows a maximum for a certain intermediate 
load while the damping coefficient steadily increases. The computed results evidence the 
paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where 
there is actually an inflow through the bearing inner diameter, accompanied by 
subambient pressures just downstream of the bearing recess edge. These results are solely 
due to centrifugal fluid inertia and advection transport effects. 

Recommendations include the extension of the computer program to handle flexure pivot 
tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft 
angular misalignments. 
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Introduction 

The importance of hybrid (combination hydrostatic and hydrodynamic) journal and thrust 
bearings and damping seal bearings as radial support elements in cryogenic 
turbomachinery has steadily grown over the past few years. Advanced primary power 
cryogenic turbopumps are compact - low count part units, operate at very high shaft 
speeds, and require of externally pressurized fluid film bearings to support the expected 
large thrust and lateral radial loads (Pelfrey, 1995). Fluid film bearings enable smaller 
and lighter turbopumps through no bearing DN life limitation and no sub-critical rotor 
operation. These mechanical elements have durability, low friction and wear, accuracy of 
positioning, and large direct stiffness and damping force coefficients. These features 
enable the design (and operation) of unshrouded impellers with a significant increase in 
the turbopump mechanical efficiency. 

The growth of an "all-fluid-film- bearing" technology for advanced and less costly (per 
launching cost) turbopumps demands the development of analytical models and design 
tools, the testing of components, and the implementation of the technology. San Andres 
(1990-1996) has perfonned the thennohydrodynamic analysis and developed computer 
programs for prediction of the static and dynamic force response of radial fluid film 
bearings for cryogenic applications. The research addresses effectively the most 
important theoretical and practical issues related to the operation and dynamic 
performance of cryogenic fluid film bearings, i.e. namely, geometric configuration, 
operating conditions, flow turbulence, fluid inertia, fluid compressibility, thennal effects, 
and two-phase flow phenomena. The computational programs are not restricted in their 
application to cryogenic fluid film bearings and seals. The codes have been validated 
with experimental data from process fluid film bearings with mineral oils, water and air 
in regimes of operation ranging from laminar flow to turbulent flows, and including the 
transition zone to fully developed turbulence. 

Further applications of compressible fluid thrust bearings in the Thennal Management 
System of the Space Station Freedom also motivate further the need for advanced models 
and computational tools. Current thrust bearing computational programs are based on the 
Reynolds lubrication of classical lubrication, lack a comprehensive thennohydrodynamic, 
do not include real fluid properties, and neglect fluid inertia (Coriolis and centrifugal) 
effects on the bearing film lands. 

Pinkus and Lund (1981) study the effects of centrifugal inertia on the perfonnance of 
hydrodynamic thrust bearings. Hashimoto (1989, 1990) extend the analysis to include 
fluid advection transport effects. These original references have been instrumental to the 
present development. Centrifugal forces, most important in the performance of parallel 
plate bearings and seals, lead to subambient pressures and may induce lubricant 
cavitation. Thus, large areas of the bearing surface could be denuded of fluid with a 
significant reduction in the bearing load capacity. Operation at high rotational speeds will 
also generate large circumferential fluid speeds with a tremendous impact on the inertial 
pressure drop at the edges of the bearing recesses. 
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A bulk-flow analysis for prediction of the static load perfonnance and dynamic force 
coefficients of angled injection, orifice-compensated hydrostatic / hydrodynamic thrust 
bearings follows. The model accounts for the bulk- flow mass, momentum and thermal 
energy transport, includes flow turbulence and fluid inertia (advection and centrifugal) 
effects on film lands and recesses, and incorporates cryogenic fluid properties using a 
NIST data base (McCarty, 1986). The computer program predicts the flow rate, load 
capacity, power loss and axial dynamic force coefficients for rigid surface, tapered land 
hybrid thrust bearings. The static and dynamic forced performance of two hydrostatic 
bearing geometries of interest for advanced cryogenic turbomachinery and oil-free 
compressors are discussed in detail. Recommendations for further research close the 
present development. 
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Analysis 

Consider the turbulent flow of a fluid within the film lands of a hybrid 
(hydrostatic/hydrodynamic) thrust fluid film bearing. The bearing geometry and coordinate 
system are depicted in Figure 1. In the following, the bulk- flow equations of motion for the flow 
on the thin film lands are detailed along with a perturbation analysis for description of the 
equilibrium flow (zeroth-order) and perturbed flow (first-order) due to shaft small amplitude 
axial motions. 


i 

H: film thickness 

t 




Figure 1. Geometry of a hydrostatic / hydrodynamic thrust bearing 

Bulk-flow equations on the film lands of a hydrostatic thrust bearing 

Flow turbulence, fluid inertia and compressibility effects are important on the thin film lands of a 
hydrostatic thrust bearing. The model assumes the fully developed turbulent bulk-flow of a fluid 
whose material properties depend on its local thermophysical state of pressure and temperature. 
The equations of mass, radial and circumferential momentum, and energy transport for the 
bulk-flow velocities, pressure and temperature on the bearing film lands are given as (San 
Andres, 1994, 1995): 

continuity: 

d(ph) | 1 8(RpHU r ) | l B(pHU e ) Q 

dt R dR R dO 

radial momentum: 

e( P HU r ) + 1 e(R P HU> r ) +l S( P HU,M e )_ l SP 

dt R dR R do R dR 
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circumferential momentum: 
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where Qbs = Qb + Qs = H B (T-T B )+Hs(T-Ts) is the convected heat flow through the bearing and 
runner (shaft) surfaces. P and T are the fluid pressure and temperature, and (U r , UQ are the bulk- 
flow velocities in the radial and circumferential directions, respectively. Refer to the 
Nomenclature for a description of all other variables. 


The wall shear stress parameters are determined as local functions of turbulent friction factors, 
and which depend on the bearing and shaft surface conditions and the flow Reynolds numbers 
relative to the rotating shaft and stationary bearing surfaces (Hirs, 1973): 
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where the turbulent flow shear factors ( k, = i< 2 , t<s) are defined in terms of the friction factors (/) 
and Reynolds numbers (Re) relative to the stationary bearing ( B ) surface and the shaft (5) 
rotating surface. The functional forms are: 


k , 


= K e =-{ k b +k S ); k B =/« Re s ; «s =fs Re s 


where 


Re B = 


r pH_] 


^U 2 r +U 2 0 ; Re s =[^\^U 2 R +(U e -nR) 


v M j 


(3.a) 


(3-b) 


The cryogenic liquid properties (A, :,C P ) are extracted from the Benedict-Web-Rubin equation of 
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state as given in the standard data base of McCarty (1986). 

Bulk-flow equations at the recess of a hydrostatic thrust bearing 

In a hydrostatic bearing, the fluid at pressure (P s ) and temperature (T s ) is supplied through orifice 
restrictors into the bearing pockets or recesses. The continuity equation at a hydrostatic recess 
establishes a balance among the mass flow through the feed orifice (Mr), the flow through the 
boundaries of the recess into the film lands (My), and the accumulation of fluid mass within the 
recess volume, V R =[A R (f I+! Ir) + V supp i v ] . The conservation of mass flows at a bearing recess is 
given as, 


Mr. - M rj + P R 


dt 


dP 


dP 


+ PrVrijPp ~ Pt ~ , r — , 


dt 


dt 


R, 


( 4 ) 


where (P, T) R are the averaged pressure and temperature within the recess, and (p, ftp ,ftr J r are the 
fluid density, fluid compressibility coefficient and thermal expansion coefficient, respectively. 
Recall that the fluid material coefficients are defined as 


ft P -+ 


1 dp 

p dP 



1 dp 

p dP 


p 


( 5 ) 


The mass flow (Mr) through the orifice is given in terms of the pressure drop from the feed line 
at supply pressure (P s ) to the recess inlet pressure (Pr), i.e. 


M Ri -C d . A 0 . 2 p Roj (p s 



-,Nrec 


( 6 ) 


7T i . . 

where A 0 =—d~ is the orifice area and C d is an empirical discharge (loss) coefficient. The mass 
flow rate leaving the /^-recess through its boundary (P) and into the thin film lands is given as 

M ri =lpHU-f) dr t (7) 

with the fluid velocities evaluated from the governing equations on the film lands. 

A global flow energy balance is also accounted for at the bearing recesses. The simple model 
includes the mechanical energy dissipated by viscous shear, the heat-carry over (advection) from 
upstream conditions and the thermal mixing effects. The transport of thermal energy equation 
within a bearing recess is written as: 

(c t v R =C r M R T,~ iC r T(pHu)ijdr+T os n ) ( 8 ) 

\ dt r ‘ 1 1-1 
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where T q r is the shear viscous drag torque on the recess area. 


The analysis also considers the circumferential bulk-flow and pressure rise within the hydrostatic 
recesses as a one-dimensional Rayleigh-step bearing. The simple model also includes the 
momentum transfer due to angled fluid injection against shaft direction. San Andres and Childs 
(1997) provide full details on this analysis. The edge recess pressure ( P R ) just before the inlet to 
the film lands is given by the approximate relationship (San Andres, 1995): 


Prc- ~ Pr + Mr 


2{H + H R )~ 


K e U e ~K s 


QR 


\ 

Jr 


( 9 ) 


where bn is the circumferential length of a hydrostatic recess 1 . 


High-speed thrust bearings can develop significant centrifugal flow effects. San Andres (1997) 
presents a simple model for prediction of the radial pressure rise within the recess region due to 
the centrifugal fluid acceleration. The radial pressure within the recess is given by: 


P-P R + P R(lyn 


/ „ \ 


( u J 

2' 

(s 2 -i) 

1 + 

>R 




V SU 0 R J 



2 P k,.U R 

> ln(s ) 


p(H + H R ) ! U, 


2 

e R 


J / i=l,...,N n 


( 10 ) 


where (s=R/Rr), and (U rR ,U 0r = a Q R r ) are the radial and circumferential bulk-flow velocities 

1 , 

within the recess, and P R = — p U e is the dynamic pressure head due to centrifugal fluid 

dyn 2 K 

inertia. The expression above shows the variation of recess pressure due to the following effects: 

a) viscous shear decreasing the pressure as the radius grows, i.e. towards the outer side of the 
bearing. 

b) centrifugal forces due to fluid rotation which raise the pressure towards the outer radius of 
the recess, and 

c) advection of fluid momentum in the radial direction that decreases the pressure as the radius 
within the recess grows. 

The viscous and centrifugal fluid inertia effects on the generation of the recess pressure field, 
equations (9) and (10), are superimposed. This assumption is grounded on simplicity rather than 
on a thorough analysis of the complex flow field within the hydrostatic pockets. 

The local acceleration of fluid from a deep recess into the thin film lands causes a sudden 
pressure drop. The pressure [p R( )at the entrance to the film lands is modeled from Bernoulli’s 
equation as, 


1 Note that this length varies in the radial direction. 
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( 11 ) 
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where E> = {<f , c„ , g (hl , g a/ } are empirical entrance loss coefficients for the edges of a recess at its 
inner radius, outer radius, upstream and downstream circumferential directions, respectively. The 
sudden pressure drop is accounted for only if the fluid flow effectively enters the thin film lands. 

Governing equations in dimensionless form 

Let dimensionless coordinates and flow variables be defined as: 
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where the (*) denotes characteristic values, and with U, = v 5 — as a characteristic 

ju , R. 

pressure flow velocity. Dimensionless orifice and recess outflow mass flow rates, recess area and 
volume are defined as 


M R,r ~ 
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The bulk-flow equations of motion on the thin film lands become in dimensionless form: 

continuity: 

j. djph) | 1 d(rp hu r ) | 1 d{p hu g ) _ Q 
dr r dr r 80 

radial momentum: 


Re p * 


a 
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(13 .b) 


circumferential momentum: 
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(13. c) 
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energy transport: 
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(13. d) 
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where /t= ~ ^ and cr = ^ ^ are the shaft speed and frequency numbers. 
C* 

— is a nominal Reynolds Number due to pressure flow, and 
R* 


Re P * = 


U, 

p* PI * c > 




M* ) 

f 


Re a = Re P * cr = 


p - ( ') •• i 




U* 

is a squeeze film Reynolds number. E c = is the Eckert heat 

T* C p* 


v P* J 

transfer number, and {h B , I! s ) are dimensionless convection heat transfer coefficients. The 

f Re p . * 

ratio — — may be interpreted as the effect of heat convection relative to shear dissipation. 

I E C J 


In dimensionless form, the mass conservation equation at the /,/, -hydrostatic recess takes the 
form, 


Mr — Mr + cr p R 


Pjl + v h s -P-b A 

dr R '\ Pp dr Pl Sr 


i=l Jf n 


(14. a) 


where M R .=8 Ro . 


pRoi [ l - p Rl T 2 > an d M r . phu-rj dT, 


(14- b) 


with an orifice parameter defined as, 



( 15 ) 
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The energy transport balance within a bearing hydrostatic recess is written as: 


C P V R, 


^ d(pT ) R ^ ^ ^ 
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The circumferential and radial pressure rises within the 6/,-hydrostatic recess are given in 
dimensionless form as: 
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while the entrance pressures to the film lands bounding the /,/,-hydrostatic recess are expressed 
as, 





- P R e + 


Re p * P, 


KPe ) 
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R J 


l r,e 
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(18) 


Film thickness expression and boundary conditions 

The thrust hydrostatic/hydrodynamic bearing maybe composed of: 

a) a single continuous (360E) pad with (N rec ) recesses distributed around the bearing area, or 

b) a number of pads (N pa( j) separated by radial grooves. Each pad may contain one or more 
recesses. 


For a bearing pad of angular extent (l pa d) and with its leading edge at the angular location 
(h_pad), the film thickness including circumferential ( *2) and radial (*r) slopes is given by 
(Pinkus and Lund, 1981): 
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( 19 ) 


(j_j \ (j_j \ 

where, referring to the geometry presented in Figure 2, 56 =- — l —^ — — ; Sr =- — ^ — — , are the 
slopes (taper) in the circumferential and radial directions, respectively. 


K-' = 
edge. 


1- 


1 ( R out ~ R in) 

2 R„. f 


with C=H 2 as the film thickness at the midpoint of the pad trailing 



Figure 2. Description of a bearing pad and tapered film thickness. 


The pressures at the inner and outer bearing radii are specified as, 
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( 20 ) 


P( r in) = P D in ;P(r out ) = P Dout 

These pressures are regarded as uniform or constant, and consequently, their variation under 
dynamic shaft motions is nil. 

In a 360E bearing, the fluid pressure, temperature and velocities are singled valued in the 
circumferential direction, i.e. 

P ,T ,u r ,u 6 1 ( r,6,t)=P ,T ,u r ,u 9 ( r,6 + 2x,t ) (21) 

while for a bearing pad, the pressures at the leading and trailing edges are given as a (known) 
function of the pressures (P D . ,P D ) specified at the inner and outer radii, respectively. A ram- 

pressure effect due to fluid inertia at the leading edge of a bearing pad is also considered as per 
the analysis of San Andres (1996). 

Note that the momentum and energy equations on the film lands are of hyperbolic character, and 
consequently, no boundary conditions are required for the discharge temperature and fluid 
velocities at the bearing inner and outer radii. The balance of flow leaving the bearing recesses 
and entering the film lands provides the (inner) boundary conditions for the velocity fields on the 
thin film flow region, see equations (18). 

Perturbation analysis of the flow field 

Determination of the bearing dynamic force coefficients requires of a perturbation analysis of the 
flow field for infinitesimally small shaft axial motions about an equilibrium position. This 
elegant procedure renders sets of zeroth- and first-order flow equations for evaluation of the 
bearing static load capacity, shear torque, bearing and recess flow rates, and the dynamic force 
coefficients. 


Consider small amplitude ()z) axial motions of the thrust collar (shaft or runner) at frequency (T) 
about an equilibrium position denoted by (ho). The film thickness is expressed in dimensionless 


form with h, = 


h = h 0 +hj e‘ 


f A 


c* 


as, 


Sh_ 

dz 


= i hj e' 


i = V-7 


( 22 ) 


All the flow variables, pressure, temperatures and velocities are also expressed as the 
superposition of zeroth-order and first-order flow fields describing the equilibrium and the 
perturbed fields, respectively, i.e., 


</> = <f>O+0i 


<?>={ 


u„ ,u, 


,P ,T , p ,/J ,K r ,K e ,K s ,etC 


.} 


(23) 


2 The inner and outer radii pressures may be different depending on the closeness of the bearing to the discharge and 
suction zones of a centrifugal impeller in a cryogenic turbopump. 
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Substitution of these definitions into the thin film land equations (13) leads to the zeroth- and 
first-order governing equations for the fluid flow: 

Zeroth-order bulk-flow equations on the film lands 

continuity: 

1 d{rp 0 h 0 u r0 ) | 1 d(p 0 h 0 u eo ) _ Q 
r dr r 56 

radial momentum: 
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circumferential momentum: 
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energy transport: 
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First-order bulk-flow equations on the film lands 

with hi=l for pure axial shaft collar motions, i.e. without shaft angulations or dynamic 
misalignments: 

continuity: 
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(24. a) 
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radial momentum: 
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(24. b) 


circumferential momentum: 


h dP ) 
r d6 

Re * — 
r 


- Yghhj + Y 9r U rl + Yog U 91 + Yop P l + Y9T T 1 + 
d(rp o h o u rO U 0 j) | d{p 0 h 0 u 91 u e0 ) | 


Po h. 


dr 
r 
o 


do 


ujr + u r0 + 


du 


90 


do 


u 91 + Po h() 


r^+u 


dr 


90 


l rl 


(24. c) 


energy transport: 
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The first-order wall shear stress coefficients ((‘s) are given in Appendix A. 

Zeroth- and first-order flow equations at a bearing recess 

Perturbation of the recess mass flow and energy transport equations proceeds with the following 
linear combination of equilibrium and dynamic fields, 
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(25) 


M R -M r 0 +M R] e ; M r —M rn +M ri e ;P R —P Rn +P Rl e ;T R —T Rn +T R ,e 
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to render the zeroth and first-order equations for mass flow conservation at each recess, i.e., 
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(28. a) 


and M r = $ r {ft, h 0 u 0 + p 0 h, u 0 + p 0 h 0 u, )•?; dr \ ; i=1 


w 


(28. b) 


are the first-order mass flow rates through the recess boundaries into the film lands. 

The first-order energy transport equation at the recess and the pressure rise/drop equations at the 
recess edges are omitted for brevity. 

Fluid film reaction force and axial force coefficients 

The fluid film force (F z ) reacting to an applied external load (IV/) is given by integration of the 
hydrodynamic/hydrostatic pressure field on the bearing surface, i.e., 


F z =-W z =UP-PJRdRd0 


(29) 


Substitution of the calculated zeroth- and first-order pressure fields into the equation above leads 
to the determination of the bearing load capacity and the dynamic stiffness (K/z), damping (C/z), 
and inertia (M Z z) axial force coefficients, i.e. 


Fz = (Ps-PJ^Bh, P o’-drd0 


(30) 


K/z 0) M zz +ico C zz - ^ ( P s P a ) Ar fa B Pi r dr dd 


(31. a) 
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In general, bearings operating with incompressible fluids lead to force coefficients independent 
of the excitation frequency (1). However, compressible fluids and bearings with recesses of large 
volume give rise to force coefficients that are complicated functions of the frequency (I). In this 
case, the results from the model should be interpreted as frequency dependent impedance 
functions whose real and imaginary parts render dynamic stiffness and damping coefficients, 
respectively, i.e., 

K zz {(o)+i(o C zz (co) = ~( P s - p a) A B k B Pi rdrdd ( 3Lb ) 

San Andres (1991) discusses the effects of recess volume - fluid compressibility on the dynamic 
force response of hydrostatic bearings operating with compressible fluids (such as LH 2 ) and 
provides criteria to avoid pneumatic hammer (loss of damping) instability. It is important to note 
that the notion of frequency dependent force coefficients does not imply these coefficients are 
nonlinear in terms of the dynamics of a rotor-bearing system, i.e. the force coefficients are not 
functions of the amplitude of motiorn . 

The shear drag torque (T 0 ) on the rotating member is generated by the shear stresses acting on 
the runner surface, and given as 

I" R 2 dRdgJt'PPL ) f \L!l^+JL{ k u K ( u Ar )\ r ! drde (32) 

a b { c j j B [2 r d6 4 h J 

Numerical method of solution of flow equations on the film lands 

A control-volume finite difference scheme is implemented to solve the coupled, nonlinear 
PDEs of mass, momentum and energy transport governing the fluid flow in a 
hybrid/hydrostatic thrust bearing. The flow field is represented by a series of discrete nodal 
bulk-flow velocities, fluid pressures and temperatures on staggered grids (Patankar, 1980). An 
advantage of the staggered grid is that the potential unrealistic wavy pressure solutions arising 
from a single grid are avoided. Algebraic difference equations are obtained by integration of 
the governing equations on finite size control volumes and using the full advection SIMPLEC 
algorithm of Van Doormaal and Raithby (1984). The effectiveness of the numerical algorithm 
has been demonstrated on the turbulent bulk-flow analysis of cryogenic fluid film radial 
bearings and seals (San Andres, 1990, Yang, 1992, Yang et ah, 1994). 

Let the dimensionless zeroth- and first-order bulk flow field variables \u r ,u 9 ,P ,t} 0 l be 

represented by the discrete fields jlJ, V, P, T } and ju^V^P^T,}, respectively. The continuity 
equation, circumferential and radial momentum transport equations, and the energy transport 
equation are integrated on the P-, U-, V-, T-control volumes depicted in Figure 3 to render sets 
of algebraic difference equations. Appendix B details the steps in the numerical procedure and 
includes all the equations derived and as implemented in the computer program. 

3 Of course, this last assertion is only valid for truly small amplitude motions about an equilibrium configuration. 
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P,T-control volume 



Figure 3. Staggered control volumes for integration of bulk-flow equations in 
polar coordinates (r,0) 


The governing difference equations on the bearing film lands are solved iteratively. An 
approximate solution for the flow field variables is used to initiate the calculations 4 . The 

{ sje I • • • 

U ,V (are first calculated from the circumferential and radial 
momentum transport equation using an estimated pressure field (p*). Next, the pressure 
correction equation based on the continuity equation provides a correction pressure (p') and 
velocity fields {u',v'}. Solution to the energy transport equation renders an estimate of the 
temperature field (t). Updated pressure and velocity fields, 
U = U +u';V = V +v'; P = P +p are substituted again into the momentum equations 
resulting in (generally) more accurate velocity and pressure field solutions. The iterative 
procedure above is repeated until the sum of (absolute value) residuals of mass flow on the 
control volumes is below a prescribed tolerance value, typically 0.1% of the global mass flow 
rate on the whole domain of solution. See Appendix B for full details on the numerical 
algorithm. 

The numerical method proceeds to balance the flow and momentum equations in the recesses 
of a hydrostatic bearing once a solution to the fluid flow equations on the film lands is 
obtained with appropriate boundary conditions. A Newton-Raphson algorithm provides 

4 The approximate analytical model corresponds to a highly idealized geometry and overly simplified flow 
condition. Nonetheless, the approximation provides a reasonable well posed field which allows rapid convergence 
to the final solution. 
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improved estimates of the recess pressures and velocity fields. The procedure is continued 
with the solution of the equations on the film lands until convergence is obtained in both film 
lands and recesses. 


The iterative scheme for balance of flow at bearing recesses. 

The recess pressures on each pocket of a hydrostatic bearing are updated after the computation 
of the flow field on the film lands. This procedure is performed iteratively until the set of mass 
flow rates into the film lands is balanced by the mass flow through each of the feed orifices in 
the recesses. Recall that the zeroth-order mass conservation at each recess is given by (the 
subscript “0“ is omitted for brevity): 



(33) 


At an intermediate step in the calculation of the flow field a set of “old” values of recess 
pressures and mass flows into the film lands is known, i.e. \P R ;M r \ (old 1 . A Newton- 

‘ 1 i=l N rec 

Raphson scheme allows the determination of “new” or improved recess pressures from the 
balance of the flow out from the orifices with the flows into the film lands. A first-order Taylor 
series expansion of the mass flow equation at each recess gives, 
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S ‘ k e p r 


M r i {p^ ,d > + AP Rk ) M P (pg*>) 

U 1 — , i,k = l,2,...,N r 


ap d 


(34) 


(35) 


as a set of mass flow (from recess to land) influence coefficients calculated with small (finite) 
pressure increments \AP R ) around the (old) values of recess pressure. The fluid density 

is also a function of the flow pressure and temperature; and a Taylor series expansion leads to, 


pOiew) -(old) d P Rj S p.j 

1 ' d P R; d T Ri 


= pf> + (3p R ® AP r . - j8%fAT Ri J, i = l, 2,...,N n 


(36) 
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where [aT r = T R ew> - T R ld> ). =/ N are variations in the recess temperatures. Substitution of 

the expansions into the recess mass flow rate equations for each recess leads to the following 
linear system of equations: 


v rec 
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S ik + 5 ik S Ro k 


l-P PRk (l-P Rk )-Pr Rk AT Rk 


(old)p (new) 
R k 


i.e, 
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where 8 


ik 


1 , if i = k 
0, otherwise 

iterative process until (m f 


, and from which updated (new) recess pressures are obtained in the 



The numerical procedure of the iterative scheme is described briefly as: 


a) Given initial recess pressures [p£° ld ' ) . 1 1 N > calculate the discrete flow fields 

[U*,V*,P*, T*, and {t r . ). ] on the thin film lands of the bearing. 

1 i=l,-,N rec 


b) Given a pressure increment to the k t h recess (A P Rk ) while all the other recess pressures 
remain unchanged, P R old) , i p k , recalculate the corresponding thin film land flow fields, 

and detennine the mass flow influence coefficients, , l . , ,, , and \ATn. , for 

’ L lki i,J=l,-,N re c ’ V R i >i=l,..,N rec 

all recesses in the bearing or pad. 


c) Reset the flow fields to the original ones, {U*, V* P * , T*, and (t r . )._j n calculated in 
step (a) and return to step (b) until all the recesses ( k = 1, 2,..,N rec ) are covered. 


d) Solve the system of equations (37) to update the recess pressures to \P R " ew) ) . 12 N • 

e) Return to step (a) and replace the initial recess pressures [P R old) ] } 2 N by the updated 

ones ( p( new ) The iterative process continues until the difference of the orifice 

V R ‘ ' i=J,2,...,N rec 
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mass flow, [M , and flow rate into the film lands, \M R .) , for each recess 

’V Ri)i=l,2,...,N r ec V 1 i=U ^rec 

is within a tolerable value (about 0.6%). 


Numerical solution of the first-order bulk-flow equations 

Once the solution to the zeroth-order equations on the bearing film lands and hydrostatic 
recesses has been obtained, the (perturbed) first-order flow fields {(J , ,V t , P t , T , } are calculated 
for a given frequency (T). 

The first-order equations are linear with coefficients fully determined by the zeroth-order flow 
field. The linearity of the equation enables a solution of the first-order flow field to be found by 
the component (superposition) method introduced by San Andres (1990). In this procedure, the 
first-order flow field solution is given by the form, 

</) (38) 

i 1 

The first part in the components field, (f) 'J , is the solution to the non-homo geneous fonn of the 
first-order equations with the first-order recess pressures, \P Rl ( =0 . The second part 

i i-l,..,N rec 

corresponds to the addition of the component solutions \(j)] }. =; N determined from solution to 
the homogeneous form of the first-order equations, i.e. with {/ 2 y } = 0 , and a complex unit value 
for each recess pressure, i.e., \p Rj j =( 1+i ). 

1 i=l--N re c 

The numerical scheme for the first-order equations on the film lands follows closely the one 
adopted for the zeroth-order equations. The same steps are followed for the solution of both the 
homogeneous and non-homogeneous first-order flow equations. Once all field components have 
been calculated, the first-order flow field is substituted into the first-order mass flow equations 
for each recess, equation (27), and a linear system of complex algebraic equations is obtained for 
calculation of the first-order recess pressures )P R , i.e: 




r « + C l TL *S m Ci )p. =-(m° +i*p A. h„ *C T „ Ti ); , 
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where 
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( 40 ) 


M r . — jy-. ( Po h 0 Uj + p 1 h 0 u 0 ) • n dr i , 

are the 0 t h- and A: t h-component first-order flow rates into the film lands at the boundary of the /th- 
ree ess (9j). 

The first-order mass conservation equation (39) at the recesses has the same form as the zeroth- 
order equation (27). However, the first-order equation represents a system of algebraic equations 
in the complex number domain, while the zeroth-order equation contains only real numbers and 

variables. Note that the set { P }k=i,...Nrec in equation (39) needs to be evaluated only once since 
the first-order equations are linear and the superposition principle applies. On the other hand, 
{Pp" ew) } k =i,... N r ec in the non-linear equations (27) are calculated iteratively until the recess flows 
balance the orifice flows. 

Once the zeroth-order and first-order flow field are obtained, the fluid film forces and the 
dynamic force coefficients are calculated by integration of the zeroth- and first-order pressure 
field over the bearing surface, see equations (30) and (31). 
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The computer program hydrothrust 

The computer program developed in Fortran 90 is named hydrothrust and runs as a 
console application in a personal computer under the MS Windows operating system. 

The software package includes a windows based help fde ( hthrust.help ) and several 
examples featuring the options and capabilities of the program. 

hydrothrust calculates the static and dynamic force performance characteristics for the 
following bearing types: 

1 . hydrostatic / hydrodynamic thrust bearings with orifice compensation, 

2. annular face seal with a pressure drop from inner diameter to outer diameter, 

3. plain hydrodynamic thrust bearings. 

hydrothrust includes the following thermal models: 

adiabatic surfaces, i.e. insulated shaft and bearing surfaces, 
isothermal shaft at specified temperature and insulated (adiabatic) bearing, 
isothermal bearing at specified temperature and insulated (adiabatic) shaft, 
isothermal shaft and bearing surfaces. 

isothermal shaft and radial heat flow through bearing (stator), 
adiabatic shaft and radial heat flow through bearing (stator). 

hydrothrust provides numerical predictions of: 

• bearing flow rate or seal leakage 

• friction torque, power dissipation and temperature rise, 

• load capacity if bearing minimum film clearance is given, or bearing film clearance if 
the external thrust load is given. 

• axial stiffness, damping and inertia force coefficients evaluated at a specified 
excitation frequency. 

with the following (single phase) fluids: 

(1) parahydrogen, (2) oxygen, (3) nitrogen, (4) methane, (5) water, (6) oil, (7) air, 

(12) baro tropic fluid. 

A detailed description of the program operation and input/output calculation options is 
given in the help file hthrust.help. 

Numerical results and discussion 

The analysis and computational program are applicable to a wide range of thrust bearing 
applications including low speeds and pressures with viscous mineral oils (laminar flow 
bearings). However, the present development constitutes an advanced model for high 
speed hydrostatic / hydrodynamic thrust bearings for implementation in modern compact 
cryogenic liquid turbopumps. Unfortunately, experimental results for this application are 
not (yet) available in the open literature. Yet, undoubtedly, the interest on this novel fluid 
film bearing technology will render shortly important advances and actual applications. 

Numerical predictions for thrust bearings operating in the laminar flow regime are given 
as examples in the hydrothrust software bundled package. Comparisons have been 
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performed with numerical results of Hashimoto (1989) and Zirkelback and San Andres 
(1996) for hydrodynamic thrust bearings. These examples are omitted here for brevity. 

Two high speed hydrostatic/hydrodynamic thrust bearing applications are hereby 
discussed at length. The first one corresponds to a thrust bearing geometry for the 
advanced liquid hydrogen ( ALH) turbopump at nominal operating conditions, and the 
second example refers to a dual application hydrostatic thrust bearing for implementation 
in commercial (oil free) R134a compressors. 

Example of a hydrostatic thrust bearing for a cryogenic turbopump 

Table 1 lists the geometry and operating conditions of a liquid hydrogen, eight recess 
hydrostatic bearing for an advanced cryogenic turbopump. The pump design is uniquely 
compact with a low number of parts, and relies on hydrostatic radial and thrust fluid film 
bearings for accurate positioning of the rotor. Sound design and reliable operation of the 
fluid film bearing elements also allow for unshrouded impellers with a significant 
increase in the pump mechanical efficiency. 

Table 1. Hydrostatic thrust bearing for an advanced cryogenic fluid turbopump 


Geometry, N rec = 8 

SI dimensions 

English dimensions 

Inner diameter, D,„ 

48.26 mm 

1.9 inches 

Outer diameter, D out 

99.06 mm 

3.9 inches 

Recess diameter, D R 

73.66 mm 

2.9 inches 

Recess radial length, L R 

10.16 mm 

0.4 inches 

Recess arc length, 1 R 

22 E 


Recess depth, H R 

0.254 mm 

0.0100 inches 

Film clearance 

0.038 mm 

0.0015 inches 

Recess/Bearing area ratio 

0.25 


Orifice diameter 

Varies 


Empirical parameters 

Orifice discharge coefficient, C d 0.90 


Entrance loss coefficients, 

0.00 


' > ri, ^ro, ' > 2u , >2d 

Inlet swirl coefficient, V 

0.50 


Reynolds numbers 

Re c =162,740 

Re®=125.2 

Operating conditions 

SI units 

English units 

Speed (Z) 

17,802 r/s 

170,000 rpm 

Surface speed, ZR,„ , XR,„„ 

442 - 908 m/s 

1450 -2979 ft/s 

Fluid: liquid hydrogen, 

Supply temperature, T, 

45 EK 

81 ER 

Supply pressure, P s 

393 bar 

5,700 psia 

Exit pressure, P =P Di „=P ' Do ut 

81.6 bar 

2,200 psia 

Fluid properties 
(supply/discharge) 

Density, [A s , A a ] 

(81.6, 65.3) kg/m 3 

(5.09, 4.08) lb/ft 3 

Viscosity, : a ] 

(0.0173, 0.0095) mPa.s 

(0.0025, 0.0014) /Reyns 

Sound speed (v s ,v a ) 

(1793, 1299) m/s 

(5882, 4262) ft/s 

Fluid bulk modulus (1/3 P ) 

1,206 bar 

17,495 psi 
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The operating conditions are extraordinary in tenns of the high rotational speed (170 
krpm) and available supply pressures. Note the large values of the shaft collar surface 
speeds (PR), ranging from 442 to 908 m/s for the inner and outer radii, respectively. 
Furthermore, the large magnitude of the circumferential flow Reynolds number, 
Re c =A s PR out C/: s = 162,740, denotes a highly turbulent flow regime with dominance of 
fluid inertia effects. 

The analysis example considers a fixed clearance (C) equal to 38:m and intends to 
determine the orifice diameter necessary to achieve the maximum axial thrust stiffness. 
The numerical computations are carried out for increasing values of the recess pressure 
ratio, p ra tio = (PR-Pa)/(Ps-Pa), spanning from 0.1 to 0.9. The optimum stiffness condition is 
achieved at a particular recess pressure ratio, typically around 0.50 for laminar flow 
conditions. The example assumes that the turbopump has two thrust bearings acting in 
parallel, both balancing the thrust from the uneven axial pressure distribution in the 
impellers. The author has not been able to find (yet) definite information regarding the 
magnitude of the thrust loads at the operating point. 

The numerical computations demonstrate the paramount effect of fluid inertia on the film 
lands of the thrust bearing as well as at the edges of the bearing recesses. The following 
figures depict computed results from three fluid flow models, 

a) full inertia , includes fluid inertia (advection and centrifugal) at the film lands and 
recess edges, 

b) recess edge inertia, accounts only for the sudden pressure drop at the recess 
boundaries as the fluid flow accelerates towards the film lands, 

c) no fluid inertia , i.e., equivalent to the solution of the turbulent flow field within the 
film bearing lands as given by the classical Reynolds equation. 

Figure 4 depicts the calculated orifice diameter versus the p rati0 for the thrust bearing, and 
determined from the balance of flows through the orifices and the flow into the film 
lands. The model without fluid inertia renders the largest orifice diameters denoting 
larger flow rates (see later Figures 7) since fluid inertia acts as an additional flow 
resistance, in particular at the recess edges. Figure 5 shows the axial fluid film bearing 
reaction force versus the recess pressure ratio. The full inertia model renders 
approximately 23% less load capacity than the inertialess fluid flow model. As the p mt io 
increases (larger orifice diameter) the load also increases proportionately. However, 
larger loads are also accompanied by larger flow rates and lower stiffness coefficients, as 
will be seen shortly. The bearing reaction force is approximately equal to 31.6 kN (7,100 
lbs) at a p ratio- 0.55 for the full inertia model. The significant load magnitude is (tough) 
23% of the bearing specific load given by (P s -P a )A B . 

Figures 6 depict the calculated pressure profiles for the thrust bearing as determined from 
the full inertial model. The top figure shows the pressure at the radial location 
corresponding to the recess (and orifice) diameter versus the circumferential coordinate 
and for decreasing pressure ratios. The bottom figures display the pressure contours and 
pressure surface over the bearing for a p ra tio = 0.50. Note the large (sharp) drops in fluid 
pressure at the recess edges (downstream and radial). Note that the pressures in the film 
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lands are quite low, thus indicating a dominance of fluid inertia effects. Furthennore, for 
the lowest pressure ratio, p ra tio = 0 . 1 , the minimum pressure within the film lands is lower 
than the discharge pressure (151.7 bar) due to the rapid fluid acceleration as it leaves the 
recess volumes. 

Figures 7 display the bearing mass flow rate (M) versus the recess pressure ratio 
(increasing orifice diameters). The bearing flow rate increases steadily as the p ratio 
increases, and the full fluid inertia model shows the smallest flow rates. Figure 7(b) 
depicts the mass flow rates leaving the bearing through the outer (M Rout ) and inner (M Rin ) 
radii, and where M=M Rout +M Rin . Note that the flow through the inner diameter is much 
lower than the flow through the outer radius. The difference is most obvious for the full 
inertia model due to the centrifugal fluid inertia effect. 

Figures 9 through 1 1 show the stiffness (K zz ), damping ( C zz ), and inertia ( M zz ) axial 
force coefficients versus the recess pressure ratio for the three flow models. Note that 
only the full fluid inertia model renders inertia force coefficients. The stiffness 
coefficient ( K zz ) has an optimum value for a p rat i o ~0.5-0.65 depending on the fluid flow 
model. Note that the full inertia model provides the lowest stiffness magnitudes while the 
recess edge inertia model predicts an optimum force coefficient for a larger recess 
pressure ratio. The damping coefficient ( C zz ) has a maximum for low values of p ra tio, and 
the full inertia flow model provides the largest damping values. The predicted results 
show similar trends as in radial hydrostatic bearings operating at high Reynolds numbers 
(San Andres, 1990). The test measurements of Franchek, et al. (1995) for water 
lubricated hydrostatic radial bearings provide further evidence on the effects of fluid 
inertia reducing the stiffness and increasing the damping coefficients for high rotational 
speed operation. 

Figure 10 shows the inertia force coefficient ( M zz ) to decrease steadily as the recess 
pressure ratio increases with a minimum value at a p rati0 ~ 0.5. The negative value is due to 
the fluid compressibility and indicates the frequency dependency of the force 
coefficients. Note that for shaft vibrations with a frequency synchronous with the 
rotational speed, i.e. at 2,833 Hz, the inertia force coefficient although small could 
introduce (in general) a marked de-stiffening effect. 

Figure 1 1 depicts the drag torque versus the recess pressure ratio for the three fluid flow 
models. The largest torque and associated power loss are detennined from the no fluid 
inertia model. Note that at a p mtio =0.50, the power loss, i.e. the product of the torque 
times the rotational speed, is equal to 97.4 kW (130.5 Hp). This large magnitude is 
essentially due to the smallness of the film clearance (0.038 mm). 

Figures 12 and 13 depict the maximum values of the bulk-flow Reynolds number and 
fluid speed within the fluid film flow region. Note that the full inertia model predicts the 
largest fluid speeds and Reynolds numbers, thus explaining the larger damping 
coefficients since the viscous dissipation of energy is the greatest. The most important 
finding is related to the very large bulk-flow fluid speeds (and Reynolds numbers) which 
are as large as the shaft speed at the bearing outer diameter (-900 m/s). Note that the 
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sonic speed for LH 2 is just 1,300 m/s at the bearing discharge planes, and consequently, 
the flow Mach numbers are around 0.50 for most pressure ratios. 

Example of a hydrostatic thrust bearing for a refrigerant compressor 

Table 2 presents the geometry and operating conditions of a six recess hydrostatic thrust 
bearings with R134a refrigerant for a commercial compressor application. The bearing 
design demonstrates the successful extension (dual use) of the cryogenic fluid film 
bearing technology to a process fluid film bearing application. In this example, the 
bearing geometry is specified including the orifice diameter 1 . The operating film 
clearances, between 12.7 pm to 101.6 pm, are calculated for a range of thrust loads at 
two operating speeds, 10 and 16 krpm, and pressure drops equal to 5.17 and 10.34 bars 
(75 and 150 psi), respectively. The predictions include the full fluid inertia and recess 
edge inertia models to demonstrate the paramount effect of land fluid inertia on the 
performance characteristics of the thrust bearing. 

Table 2. Hydrostatic thrust bearing for R134a compressor application 


Geometry, N rec = 6 

SI dimensions 

English 

dimensions 

Inner diameter, D in 

89.13 mm 

3.51 inches 

Outer diameter, D ou , 

126.8 mm 

4.99 inches 

Recess diameter, D R 

108.6 mm 

4.28 inches 

Recess radial length, L R 

11.68 mm 

0.46 inches 

Recess arc length, 1 R 

24 E 


Recess depth, H R 

0.508 mm 

0.0200 inches 

Film clearance 

0.012-0.101 mm 

0.0005-0.004 inches 

Recess/Bearing area ratio 

0.25 


Orifice diameter 

1.70 mm 

0.067 inches 

Empirical parameters 
Orifice discharge coefficient, C d 

0.80 


Entrance loss coefficients. 

0.0, 0.0, 0.0, -0.5 


->2u 5 ^2d 

Inlet swirl coefficient, V 

0.50 


Operating conditions 

SI units 

English units 

Speed 

1047.2 - 1,675 r/s 

10,000 - 16,000 1 pm 

Fluid: R134a refrigerant 
Supply temperature, T, 

311 EK 

560 ER (100 EF) 

Supply pressure, P, 

18.96 -24.10 bar 

275 - 350 psia 

Exit pressure, P a =P D i„=PDo, ll 

13.80 bar 

200 psia 

Saturation pressure. 

9.63 bar 

139 psia 

Fluid properties 
Density, A a 

1210 kg/m 3 

75.54 lb/ft 3 

Viscosity, : a 

0.000198 Pa.s 

0.0288 /Reyns 

Fluid bulk modulus (1/3 P ) 

1,820 bar 

26,667 psi 


Figure 14 shows the equilibrium film clearance (C) versus the thrust load for the two 
operating conditions at 10 krpm and 16 krpm. The full fluid inertia model predicts 
smaller film thickness than the recess edge fluid inertia model, and the bearing with the 


1 The orifice discharge coefficient selected is based on prior experiences with hydrostatic radial bearings. 
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larger pressure drop, i.e. P s =24.1 bars, shows a larger load capacity. In the following, the 
calculated results are presented in dimensionless form in an attempt to identify unique 
trends independent of the individual physical parameters. 

Figure 15 depicts the equilibrium recess pressure rati o, p ra tio = (Pr-PJ/( P s-Pj, versus the 
dimensionless load capacity, Wz=Wz/A B (P s -Pa), for the two operating conditions. The 
recess pressure ratio increases as the thrust load increases (decreasing film clearance), 
and the fluid inertia model predicts larger recess pressures since it acts as an additional 
flow resistance at the film lands. Figure 16a shows the bearing mass flow rate ( M) to 
decrease steadily as the load Wz increases (smaller film thickness and larger p ra tw )• The 
mass flow rate is made dimensionless with respect to N rec CdA 0 (V 2 p s [P s -P a ]) 1/2 
representing the orifice flows from supply to ambient conditions. The full fluid inertia 
model shows slightly smaller flow rates. Figure 16(b) depicts the mass flow rates leaving 
the bearing through the outer and inner radii, with M= M Rout + M Rin . The recess edge 
inertia model predicts an inner radius flow slightly smaller than the outer radius flow 
since the flow area is smaller. However, the full inertia model which includes the effects 
of centrifugal fluid inertia shows an outer radius flow much larger than the inner radius 
flow, in particular at low loads (small recess pressures and large clearance). In fact, the 
flow through the inner radius could be negative, i.e. the fluid is actually drag into the 
bearing due to a local suction (less than ambient) pressure at the film lands. The opposite 
effect occurs at the bearing outer diameter where more fluid is thrown away due to the 
large centrifugal inertia force. This force is proportional to the circumferential flow 
Reynolds number, Re c =(pR 0 utC/ju), which increases as the clearance also increases (see 
Figure 22 later for the magnitudes of the maximum flow Reynolds numbers in the 
bearing). The calculated results thus indicate that operation with large clearances may 
cause fluid starvation (not a full fluid film) in certain regions of the film lands. 

Figure 17 shows the dimensionless stiffness coefficient, K ZZ =K ZZ C/A B (P S -P a ), versus the 
dimensionless load W z at the two operating speeds and pressure supplies. The 
magnitudes of the stiffness coefficients are similar for both flow models, although the 
optimum stiffness occurs at a lower load ( p ra tio~0.6 ) for the full inertia model. Figure 18 
shows the dimensionless damping coefficient, C ZZ =C ZZ CQJA B (P s -P a ) to increase steadily 
for small to moderate values of the dimensionless load W z . However, a large increase in 
damping is apparent at the largest loads and due essentially to the smallness of the film 
clearance (see Figure 14). Note that the full fluid inertia model predicts much larger 
damping coefficients for the entire range of thrust loads. 

Figure 19 shows the dimensionless inertia coefficient, M ZZ =M ZZ C{2 2 /A B (P S -P a ) decreases 
as the dimensionless load W z increases. The physical values of the inertia coefficient 
range between 1.8 to 5.2 kg and could affect substantially the dynamic stiffness of the 
bearing, K zzc j=K zz -M zz PI as shown in Figure 20. A comparison with the results given in 
Figure 17 shows that the bearing dynamic stiffness K ZZ( j is approximately 20% lower 
than the value predicted by the recess edge inertia model. This reduction is a direct effect 
of the fluid inertia (centrifugal and advection) at the film lands. 
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Figure 21 shows the torque To** versus the dimensionless load W z . The drag torque is 
made dimensionless with respect to the reference formulae 


T oT = \ 1 + 0.024 


0.4 Re,. 


\(>.65'}/ J s f2 A b R 

C 


out 


obtained using the turbulent shear friction 


factor from Constantinescu (1962). Note that a reference drag torque value based on 
laminar flow conditions would render dimensionless torque values one to two orders of 
magnitude larger than the ones depicted in the figure. The calculated results show a 
similar drag torque for the two operating conditions and decreasing as the dimensionless 
load increases. Of course, the physical magnitude of the drag torque increases as the load 
magnitude increases since the film clearance is smaller. 


The maximum bulk-flow Reynolds numbers and half the magnitude of the nominal 
circumferential Reynolds number, A Re c =(pR 0Ut C/2 p ) , are shown versus the recess 
pressure ratio in Figure 22. The bulk- flow Reynolds number is based on the maximum 
calculated fluid speed in the bearing flow region. Larger Reynolds numbers are a direct 
result of the shaft speed and larger film clearance as the recess pressure ratio decreases 
(smaller loads). The calculated results show that the flow-Reynolds numbers are very 
similar to the nominal circumferential Reynolds number, and thus demonstrate the 
dominance of hydrodynamic effects. The full inertia model predicts larger Reynolds 
numbers, i.e. larger fluid speeds, due to the centrifugal inertial effects and the rapid 
acceleration as the fluid leaves a recess a flows into the film lands. 


Figure 23 shows the predicted maximum and minimum dimensionless pressures, 
(Pmax,min-Pa)/(Ps-Pa), in the bearing flow region versus the recess pressure ratio, 
p ratio = (P. r-P a)/ ( P s~P a) ■ The predicted results are of importance since they reveal 
significant features about the full fluid inertia model. First, the maximum film pressures 
are larger than the recess pressure at the orifice discharge plane since due to viscous and 
centrifugal fluid inertia effects there is a marked raise in the pressure at the recess corner 
facing the downstream edge and outer diameter. Note that for a p ra tio = l the peak film 
pressure is 40% above the supply value. Second, the minimum film pressures are lower 
than the discharge pressure ( P=0) and the lowest magnitudes are achieved at the lowest 
recess pressure ratios. The appearance of the subambient film pressures, nearly 40% 
below discharge value (recall P a = 12.8 bar, 200 psia), is a direct result of the fluid (inertia) 
acceleration just downstream of a recess edge. The recess edge inertia model does not 
predict any subambient pressures. 

At the operating conditions, the refrigerant saturation pressure is 9.63 bar (139 psia), and 
in dimensionless fonn (P sa turation-Pa)/(Ps-P a ) equals -0.81 and -0.405 for the cases with 
pressure drops (P s -P a )= 5.17 and 10.34 bars (75 and 150 psi) at 10 and 16 krpm, 
respectively. Thus, the thrust bearing operating at the largest rotational speed is likely to 
show fluid vaporization at the downstream edge of a hydrostatic recess! The current 
model is not (yet) able to model properly two-phase flow conditions, although Arauz and 
San Andres (1998) have advanced some relevant analytical work. Nonetheless, the 
predicted results caution the operation of the thrust bearing under such stringent fluid 
flow conditions. 
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Figure 24 displays the calculated pressure fields from the full fluid inertia model. The top 
plots show the pressure field at the pocket (orifice) diameter versus the angular 
coordinate for various recess pressure ratios. Note the sharp pressure drops at the 
downstream side of a bearing recess. The bottom figure depicts the pressure surface and 
contours for a p ra tio = 0. 5 at 16 krpm. Note the pressure raise towards the downstream edge 
of a recess, and immediately the large pressure drop causing a rapid acceleration of the 
fluid. 

As a note aside the present discussion, it is worth mentioning that the procedure 
implemented to render dimensionless bearing perfonnance variables (load, flow, 
stiffness, etc) gives meaningful results which could aid to a fast (and reliable) design of 
turbulent flow hydrostatic/hydrodynamic thrust bearings. 

Conclusions 

A bulk-flow analysis and computer program for prediction of the static load perfonnance 
and dynamic force coefficients of angled injection, orifice-compensated hydrostatic / 
hydrodynamic thrust bearings have been completed. Advanced cryogenic fluid 
turbopumps are very compact, operate at extremely high shaft speeds, and require of 
hybrid (hydrostatic / hydrodynamic) radial and thrust fluid film bearings for accurate 
rotor positioning. Sound design and reliable operation of the fluid film bearing elements 
also allows for unshrouded impellers with a significant increase in the turbopump 
mechanical efficiency. 

The analysis accounts for the bulk-flow mass, momentum and thennal energy transport, 
includes flow turbulence and fluid inertia (advection and centrifugal) effects on film 
lands and recesses, and incorporates cryogenic fluid properties using a NIST data base. 
The computer program predicts the flow rate, load capacity, power loss and axial 
dynamic force coefficients for rigid surface, tapered land hybrid thrust bearings. 

Two high-speed hydrostatic/hydrodynamic thrust bearing applications are analyzed. The 
first one corresponds to a thrust bearing geometry for the advanced liquid hydrogen 
(. ALH) turbopump at nominal operating conditions, and the second example refers to a 
dual application hydrostatic thrust bearing for implementation in commercial (oil free) 
R134a compressors. 

The computed predictions reveal that the ALH bearing load capacity and flow rate 
increase with the recess pressure (i.e. increasing orifice diameters) at the nominal 
operating condition. The bearing axial stiffness has a maximum for a recess pressure ratio 
~0.55, while the axial damping coefficient decreases as the recess pressure ratio 
increases. The computer results from three flow models are compared. These models are 
a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess 
edges and film lands. The full inertia model predicts the lowest flow rate, axial load 
capacity and stiffness coefficient, yet it renders the largest damping and inertia force 
coefficients. Some fluid compressibility effects are apparent in the application. The most 
important findings are related to the reduction of the outflow through the inner radius and 
the appearance of subambient pressures on the bearing film lands. 
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The performance of the refrigerant hybrid thrust bearing is evaluated at two operating 
speeds and pressure drops. The computed results are presented in dimensionless form to 
evidence consistent trends in the bearing perfonnance characteristics. As the applied axial 
load increases, the bearing film thickness and flow rate decrease while the recess 
pressures increase. The axial stiffness coefficient shows a maximum for a certain 
intermediate load (recess pressure ratio) while the damping coefficient steadily increases. 
The computed results evidence the paramount of centrifugal fluid inertia at low recess 
pressures (i.e. low loads), and where there is actually an inflow through the bearing inner 
diameter along with subambient pressures just downstream of the bearing recess edges. 
These results are solely due to centrifugal fluid inertia and advection transport effects. 

Recommendations 

Further research is recommended to extend the analysis and continue the software 
development to include bearing geometries and operating conditions more representative 
of actual practice. In addition, the enhanced computer program should aim to satisfy the 
needs of commercial turbomachinery applications and to render a design tool of a broad 
engineering use. The flow model is currently restricted to axial shaft (collar) motions 
without angular misalignment. Thus, the analysis of the bulk-flow for shaft dynamic tilts 
in two planes is of importance to determine restoring torques and moment force 
coefficients. 

Many hydrodynamic thrust bearings are composed of rigid pads with machined tapered 
lands. However, this low-cost bearing configuration can only be optimized for a narrow 
range of operating load and (low) speed conditions. High performance turbomachinery 
requires of tilting pad hydrodynamic thrust bearings able to adapt and configure to wide 
ranges of loads and speed conditions. In these bearings the pads rotate about spherical 
pivots to form the (wedge) hydrodynamic film able to support an axial load and without 
the need of a machined tapered surface. Some particular tilting pad thrust bearings are 
also supported on complex structural (springs) elements allowing the control of axial 
angular rotor misalignments. The major disadvantages of tilting pad bearings are their 
cost, complexity in manufacturing and stack up of tolerances on assembly. Furthermore, 
conventional tilting pad bearings are not easily modified to allow for a hydrostatic feed 
and external pressurization. 

Submerged water pumps integrate a particular type of hydrodynamic thrust bearing 
where each sector pad is attached to a bearing support with a structural (elastic) radial 
web designed to provide a low stiffness restraint to pad angular motions, thus forming the 
hydrodynamic film wedge which generates the load support. Furthermore, the flexural 
webs also offer a feasible path for a hydrostatic fluid feeding port, thus allowing for an 
increased bearing load capacity and immediate application in a cryogenic fluid 
environment. The structural pad supports could also be engineered for the bearing to 
control dynamic angular rotor misalignments. 

The fluid mechanics of highly compressible fluid media operating at near sonic 
conditions and with two-phase flow generation and transition from liquid to gaseous 
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states are also regarded as important in the near future. The analyses of thrust foil 
bearings and spiral grooved thrust bearings are also noted as important for future work. 
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HTB axial load (N) Orifice diameter (mm) 


LH2 HTB: design operation, 175 krpm, Ps = 293 bar, Pi=Po = 152 bar, T=45 K, C = 38 um 



Figure 4. Orifice diameter vs. recess pressure ratio for ALH thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po = 152 bar, T=45 K, C = 38 um 



Figure 5. Fluid film force vs. recess pressure ratio for ALH thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 
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8 recess LH2 bearing: 170krpm, Ps=393 bar, Pa=151.7 bar 



(a) 


LH2 thrust bearing: 170 krpm, Ps=393 bar, Pa=151.7 bar 



R/Rout 


LH2 thrust bearing: 170 krpm, Ps=393 bar, Pa=151.7 bar 


"plotpres50" 



Figure 6. Pressure distributions for LH2 thrust bearing, (a) at pocket diameter for varying 
pressure ratios, (b) pressure contours and (c) pressure surfaces for pratio=0.50 
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Flow rate through outer and inner diameters (kg/s) Flow rate (kg/s) 


LH2 HTB: design operation, 175 krpm, Ps = 293 bar, Pi = Po = 152 bar, T = 45 K, C=38 um 



Figure 7(a). Flow rate vs. recess pressure ratio for ALH thrust bearing 

(Effects of fluid inertia at film lands and recess edges) 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po=152 bar, T=45 K, C=38 um 



Figure 7(b). Flow rates through inner and outer radii vs. recess pressure ratio for ALH 
thrust bearing. (Effects of fluid inertia at film lands and recess edges) 
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LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po=152 bar, T=45 K, C=38 um 



Figure 8. Axial stiffness coefficient (Kzz) vs. recess pressure ratio for ALH thrust bearing 

(Effects of fluid inertia at film lands and recess edges) 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po = 152 bar, T=45 K, C=38 um 



Figure 9. Axial damping coefficient (Czz) vs. recess pressure ratio for ALH thrust bearing. 


36 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po = 152 bar, T=45 K, C=38 um 



(Effects of fluid inertia at film lands and recess edges) 

Figure 10. Axial inertia coefficient (Mzz) vs. recess pressure ratio for ALH thrust bearing 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po = 152 bar, T=45 K, C = 38 um 



Figure 11. Drag torque vs. recess pressure ratio for ALH thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 
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Minimum (bulk-flow) Reynolds number Maximum (bulk-flow) Reynolds number 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po = 152 bar, T=45 K, C = 38 um 



LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po=152 bar, T=45 K, C=38 um 



Recess pressure ratio 


Figure 12. Maximum and minimum flow Reynolds numbers vs. recess pressure ratio for 

ALH thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 
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Maximum fluid speed within bearing 


LH2 HTB: design operation, 175 krpm, Ps=293 bar, Pi=Po=152 bar, T=45 K, C=38 um 



Recess pressure ratio 


Figure 13. Maximum fluid speed flow Reynolds numbers vs. recess pressure ratio for 

ALH thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 
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Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



10 krpm, 

Ps=19 bar (full inertia) 

- - O - - 

(edge inertia) 

— — 1 6 krpm, 

Ps=24.1 bar (full inertia) 

❖ 

(edge inertia) 


1200 1500 1800 

Thrust Load (N) 


Figure 14. Film clearance (C) versus thrust load (Wz) for refrigerant thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 


Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



Figure 15. Recess pressure ratio versus dimensionless thrust load ( Wz) for refrigerant 

thrust bearing. 

(Effects of fluid inertia at film lands and recess edges) 
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Dual Use 6-pocket HTB for R134a compressor, Ts = 38 C 


E 


o 


E 

o 



Dimensionless load 


Figure 16(a). Dimensionless bearing mass flow rate versus thrust load ( Wz) for refrigerant 

thrust bearing. 


Dual Use 6-pocket HTB for R134a compressor, Ts = 38 C 
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Figure 16(b). Mass flow rates (dim.) through inner and outer radii vs. versus thrust load 

( Wz) for refrigerant thrust bearing. 
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Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 


o 


E 

5 



Dimensionless load, [W z /(Ps-Pa)AB] 


Figure 17. Dimensionless axial stiffness ( Kzz) coefficient vs. versus load ( Wz) for 

refrigerant thrust bearing. 


Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



Figure 18. Dimensionless axial damping ( Czz) coefficient versus load ( Wz) for refrigerant 

thrust bearing. 
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Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Dimensionless load, [Wz/(Ps-Pa)AB] 


Figure 19. Dimensionless axial inertia ( Mzz) coefficient versus load ( Wz)for refrigerant 

thrust bearing. 


Dual Use 6-pocket HTB for R134a compressor, Ts = 38 C 



Figure 20. Dynamic stiffness coefficient ( K zzd = K zz - M zz fi 2 ) versus load ( Wz)for 

refrigerant thrust bearing. 
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Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



o 


0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Dimensionless load, [Wz/(Ps-Pa)AB] 


Figure 21. Dimensionless drag torque( To«). versus load ( Wz) for refrigerant thrust 

bearing. 


Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



0.4 0.6 0.8 

Recess pressure ratio, (PR-Pa)/(Ps-Pa) 


Figure 22. Maximum Reynolds numbers versus recess pressure ratio for refrigerant thrust 

bearing. 
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Dual Use 6-pocket HTB for R134a compressor, Ts=38 C 



Figure 23. Maximum and minimum film (dimensionless) pressures versus recess pressure 

ratio for refrigerant thrust bearing. 
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(a) 


§ 


6 recess R134 HTB: 16 krpm, Ps=24.1 bar, Pa=13.8 bar, Psat=9.63 bar 



(b) 


6 recess R134 HTB: 16 krpm, Ps=24.1 bar, Pa=13.8 bar, pratio=0.58 

Largest film pressure 


\ 


"plotPr.58" 

22 

20.5 
19 

17.5 
16 



Figure 24. Film Pressures from full inertia model, (a) pressure at recess diameter for 
Pratio = 0-93, 0.58, 0.22, (b) pressure surface and contours for p ratio =0.58 
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Appendix A 

First-order wall shear stress coefficients 

The perturbed wall shear stress coefficients ((‘s) for the first-order equations of radial and 
circumferential momentum and energy transport are given by the following expressions. 
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Qbso - Qbo + Qs 


and the following coefficients for evaluation of the Moody’s friction factors, 
a M = 0.001375 ; b M =5-10 5 ; c M =1-10 4 ; e M = j 
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Appendix B 


Revised January 7, 2005 

Algebraic bulk-flow equations on the thin film lands of a thrust bearing 


The flow field is represented by a series of discrete nodal bulk-flow velocities, fluid pressures and 
temperatures on staggered grids (Patankar, 1980). The dimensionless zeroth- and first-order bulk 
flow field variables \u r , u o ,P,t} 0 ; are represented by the discrete (algebraic) fields 

{u , V , P, t}, and ju , , V, , P, , T, } , respectively. Figure B. 1 shows the arrangement of staggered 
control-volumes for discretization of the thin film flow equations in a polar coordinate system (r, 6). 

The continuity equation, circumferential and radial momentum transport equations, and the energy 
transport equation are integrated on the P-, U-, V-, T-control volumes depicted in Figure B. 1, in 
order to obtain a set of algebraic difference equations for the discrete flow variables. 


Additional Nomenclature 

A coefficients of discrete algebraic equations 

D coefficients relating velocity to differences of pressure corrections 

h dimensionless film thickness 

m mass flow rates across faces of control-volumes 

S source terms in algebraic equations 

(r,9) radial and angular polar coordinates 

{u , V , P , T } discrete zeroth-order circumferential and radial velocities, pressure and temperature 
{u,v,p,t}j discrete first-order circumferential and radial velocities, pressure and temperature 
{ifi , V , p j correction velocities and pressures 

5r, r59 sizes of control volume in radial and circumferential directions 


Superscripts 

u 

V 

p 

T 

* 


refers to U-control volume (circumferential direction) 
refers to V-control volume (radial direction) 

refers to P-control volume (continuity and pressure correction equations) 
refers to T-control volume (energy transport equation) 
estimated flow variable, 
correction flow variable 


Subscripts 

p refers to nodal value of variable at center of control volume 

E,W, refers to values of variables at East, West, South and North nodes relative to p-value 
S,N 

e,w, refers to values of variables on east, west, south and north (face) sides of control volume. 

s,n 

1 first-order (perturbed) flow variable 


Operators 

|| o, a || = max(0,a) ; (a,b) 1 = a x b 0 +a 0 bj 
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Figure B.l. Staggered control volumes for integration of bulk-flow equations 
in polar coordinates (r,0) 
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Continuity (mass conservation) algebraic equations: 


zeroth-order equation: 

P P P p 

m e -m w +m n -m s =0 


first-order equation: 

mf e -mf w + mf n -mf s + ia(ph) P r p Sr/ S0p =0 


where 

=(phU)^5r p p ; m p = (phV) P 8r p p ; 
m s P =(p hV)f r p 80 p ; m p =(p hvKr n p 80 p ; 



(ph V) L = Pl h 0 V 0 + p 0 hj V 0 + p 0 h 0 Vj ; 
and the following notation (p h U) = pj h 0 U 0 + p 0 h { U 0 + p 0 h 0 Up 

(ph ) 1 = P 1 h 0 + P 0 h 1 



(B.l) 

(B.2) 

(B.3) 

(B.4) 
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Figure B.2. P-control volume for integration of continuity equation and 
pressure correction equation 
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Circumferential momentum transport algebraic equation: 


zeroth-order equation: 


AJ,u p = a"u e +a;u w + A“U s + A“U n + S^ p 

first-order equation: 

AjpUjp = A e U 1e + A W U 1W + A S U 1S + A n U 1n + 
where 

A 11 — A u _i_ A u _i_A U + A U +S U 

A p -rt -E + °p 

Aip = A e + A w + A s + A n + {y ee + ip hRe 0 } p r p 
Re ft (m e - m w ) u + Re p>f (phv)“ 6r p u S0£ 

with the zeroth- and first-order source terms, 


_ Se r +hp(P p -P E ) 8r p 

(B.5) 

Srs P +^(P lp -P 1E ) u 8r; 

(B.6) 


(B.7) 

8r“ 80p + 



(B.8) 
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Figure B.3. U-control volume for integration of circumferential bulk-flow 
momentum eauation 


B.4 



r p u 5r p “S0“; 


(B-9) 
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^isp — {y eh hi+Yep Pi +YeT 


r“ 8r“ 80“ + 


Vj p lYe r r p 8r p u 80“ +Re Pt (p h f p 80“ [(U n - U s )“ r p “ +U p 8r; 


and 

Yep = Ye P (8p/8P) 0 +y ep (8p/8T) 0 ; Y 0 t = Ye P (8p/8P) 0 +y 0p (8p/8T) o 
with the fluid flow advection coefficients, 


^=ReJ|0,-< 
As = Re„ || 0 ,+m“ 


m:=(phU):Sr ;; 
m‘‘ =(p h V)“ r"S d u ; 


a; = Re,. 


I! 0 ,+< 


m 


:=(phui,sr 


A u n = Re, J| 0 -m u n 


m" n ={p hV) u n r“SO u - 


where the operator || 0, a || = max(0,a) 


(B-10) 


(B.ll) 


(B-12) 
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Radial momentum transport algebraic equation : 


zeroth-order equation: 


A;v p = A^V e + A v w V w +A^V s +A;V n +s; e + h p v (P s -P p r rW 


p/ P““P 


first-order equation: 


AipVip = A^V 1E + A^V 1W + A^V 1S +A;V 1N +sr sp +h^(P ls -P lp ) v r p v 80 p v 


where 

Ap = A e + A^, + Ag + A^j + S p 


Aip =A e +A; + Ag +A; +{y n +iphRe c } p r p v 8r p v +Re p * (m n -m s ) 


P* 


with the zeroth- and first-order source terms, 


(B.13) 


(B.14) 


(B15) 

(B.16) 



Figure B.4. V-control volume for integration of radial bulk-flow momentum 
equation 
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s v - 
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VO 


r p v 5r p v 50 p ; 


Sre = Re p * (p hU 2 )p 8 r p v 50 p ; 


(B.17) 


§lSp — {y rh h l+YrP f*l + YrT } r p ^ r p $0 p + 


Uip (Yre r p V 5r p v 50p +Re ft (p h ); 5r p v [(V e - V w ) v -2U p 50^ 


(B-18) 


and 


Yrp =Yrp(dp/3P) 0 +Y rp (dp/dT) 0 ; Y rT = y rp (a P /ap) 0 +y rp (dp/dT) ( 


(B.19) 


with the fluid flow advection coefficients, 


A v e = Re„ ||0,-<||, 

A' s = Re p , || 0 ,+m] ||, 


in 


ni 


= (fihUy.gr’; Mi = Re || 0,+< ||, < = (phU)'JC, 


= (phV)] r ;Se;-,Al = Re,. II 0 -ml II, 


_Qp.20) 
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Energy transport algebraic equation : 


zeroth-order equation: 


A p T p - A e T e + A W T W + A s T s + A N T N + S p 


first-order equation: 


Aip^ip - A E T 1E + A W T 1W + A S T 1S + A n T 1n +S 


T 

IP 


where 

Ap = +j4l w + ^l +^n +Spi+ || 0,S p2 || 

A?p =Al+Al +A l +Al +{y TT +ipC„h(Re„/E c )}'r p v 8r p v 60;; 
with the zeroth- and first-order source terms, 


(B.21) 


(B.22) 


(B.23) 

(B.24) 



Figure B.5. T-control volume for integration of energy transport (fluid 
temnerature) eauation 
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(B.25) 
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(B.27) 

(B.28) 


(B.29) 


Pressure correction equations 

Let the flow variables be given as 

U=U*+u'; V=V* + v'; P=P*+p'; (B.30) 


\J l =\J t 1 +u[-,V l =V* + V 1 ;P 1 =V*+p[ 


where the (*) fields are assumed to satisfy the momentum equations but not the continuity equation. 
The (‘) fields are corrections to be detennined from satisfaction of the mass flow continuity 
equation(s), zeroth and first-order. The fields above are substituted into the momentum transport 
equations, and implementing the SIMPLEC method (Van Doormal and Raithby, 1984), the following 
relationships between the correction velocities and pressure are obtained: 


u' p =d;(p;-p^) u ; v' p =d;(p^-p;) v 

u; P =d“ (p; p - p; e ) u ; v| p =D;(pis-pIp) v 


(B.31) 
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where 


D 


U 

P 



a: 


~T.K> 


d: 



a: 


-ZA^ b 


2 A nb — A e +A w + A n +A s 

nb 


(B.32) 


Substitution of the assumed flow fields into the continuity equation leads to the following pressure 
correction equations: 

zeroth-order equation: 


ApPp - A E p E + A w p w + A s p s +A N p N +S 


first-order equation: 


ApPip - A E p 1E + A w p lw + A s p is + A N p 1N +S 


ip 


where 


A P -YA P -A p +A p +A P + A P 

nb 

Ae =(phf 6r„ p DJ ; A p =(ph) p r, p S0j d; 
Aw “(ph) p 5r p Dw ; A p =(ph)„ p r p 80 p Dj 

and the mass flow source terms are: 

p /* * * *\P 

S p=-l m e- m w+ m n- m s ) 


(B.33) 


(B.34) 


(B.35) 


(B.36) 


S 


p _ 
ip - 


i m le _m lw +m in _m i s 


) p + ia(ph)| > 


r p 8r p 80 


(B.37) 


are evaluated with the (*) fields. 

The numerical solution procedure 

The governing difference equations on the bearing film lands are solved iteratively. The steps in the 
numerical solution are summarized as follows: 

a) Initial guesses for the recess pressures and the film fields for pressure ( P *), velocities (U*, V*) 
and temperature ( T ) are calculated. 

b) The coefficients for the circumferential momentum transport equation (B.5) along the first row of 
U-control volumes (off the bearing recess middle plane) in the film lands between recesses are 
calculated, and the transport equation is solved using the tridiagonal matrix algorithm (TDMA). 
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c) 


d) 


e) 

f) 

g) 


h) 


i) 

j) 


The treatment under step b) is repeated for the radial momentum transport equation along the 
second row of V-control volumes, equation (B.13). Note that the line-by-line solution procedure 
sweeps in the radial direction from the recesses towards the bearing outer (and inner) planes. 

The coefficients of the pressure correction equation (B.33) on the first row (coinciding with the 
first U-control volume row but between the first and the second V-control volume rows) are 
calculated and the pressure correction equation is solved by the TDMA procedure. 

The velocities and pressures upstream and along the current row are corrected by adding the 
correction fields to the estimated flow fields (U*, V*, P *). 

The coefficients of the energy equation (B.21) are calculated, and the equation is solved using the 
TDMA to render the temperature field in the T-control volume row. The fluid properties, functions 
of the local pressure and temperature, are immediately updated at this stage. 

Steps b) to f) are repeated for all subsequent rows until the bearing exit boundaries (inner and 
outer) are reached. The inlet velocities at the recess boundaries and entering the film lands are 
obtained by satisfying the local flow continuity constraint at the control volumes facing the recess 
boundaries. 

The iterative procedure above is repeated until the sum of (absolute value) residuals of mass flows 


on the control volumes 


I s p p , 

domain 


I sf p 

domain 


is below a prescribed tolerance value, typically 0.1% 


of the global mass flow rate. 

Once a solution to the flow field in the film lands has been obtained, the recess temperatures are 
calculated through the global energy balance equation at the bearing recesses. 

The zeroth-order flow rates from recess into the film lands are evaluated and compared with the 
orifice flow rates. In general, a difference exists between these two flow rates, and a Newton- 
Raphson scheme is implemented to update the recess pressures and to satisfy the flow continuity 
requirement at each recess. A thorough discussion on this step is given in the main text of the 
report. 


Steps (b) through (j) are successively repeated until the difference of the orifice inlet and exit flow 
rate for each recess is below a tolerable value (about 0.6%). 


The governing difference equations on the bearing film lands are solved iteratively. An approximate 
solution for the flow field variables is used to initiate the calculations. The velocity components are 
first calculated from the circumferential and radial momentum transport equation using an estimated 
pressure field, and then a pressure correction equation based on the continuity equation provides 
corrected pressure and velocity fields. Solution to the energy transport equation also renders an 
estimate of the temperature field. The updated pressure and velocity fields are substituted again into 
the momentum and pressure correction equations resulting in (generally) more accurate velocity and 
pressure field solutions. 
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